Stability analysis of dual solutions for mixed convection and thermal radiation with hybrid nanofluid flow past shrinking/stretching curved surface

This study portrays the stability analysis and dual solutions of mixed convection and thermal radiation of hybrid nano-fluid flow past stretching/shrinking a curved surface in the presence of injection/suction conditions. A hybrid nano-fluid, in which water is used as the base fluid, copper and alumina are used as nano-particles, and the magnetic field is taken into account. The present study’s findings will provide fruitful implications for future research in the field of fluid dynamics. The bvp4c method using Matlab software is implemented to get the numerical solution of the nonlinear partial differential equation transformed into the ordinary differential equation. The behavior of the first and second solutions under the governing parameters on the curved surface of dimensionless velocity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g'(\xi )$$\end{document}g′(ξ), shear stress profile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g''(\xi )$$\end{document}g′′(ξ), temperature profile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta (\xi )$$\end{document}θ(ξ), skin friction coefficient Cfs, and local Nusselt’s number Nus were visualized in figurative and tabular form. From this the following are investigated: as the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi j$$\end{document}ϕj increase, the velocity profile for the second solution decreases, and the opposite trend is observed for the first solution. For the values of K and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda 1$$\end{document}λ1, the shear stress profile increased for the first solution, and the opposite trend was observed for the second solution, though after some interval points, the inverse of this statement was observed.For the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S,Pr,Rd,M, \;\text{and}\; \lambda$$\end{document}S,Pr,Rd,M,andλ, the upwind thermal boundary layer of the first solution is larger than the second solution.For the value of M uphill, the estimation of the absolute value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda ci$$\end{document}λci increases for both the skin friction coefficient and the local Nusselt number. In the second solution, increasing the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β, Pr, S, Ec, and M has a similar effect on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g''(0)$$\end{document}g′′(0), Cfs, -\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta '(0)$$\end{document}θ′(0), and Nus. In the first solution, increasing the values of Ec, S, and Pr on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g''(0)$$\end{document}g′′(0) and Cfs results in a decrease. The first solution has a positive eigenvalue, whereas the second solution has a negative eigenvalue. Agreement between the present analysis and literature is acceptable.

Nowadays, many researchers are paying attention to the examination of hybrid nano-fluid flow past a shrinking/ stretching surface due to several implications in both industrial and manufacturing processes.These include polymer and wire drawing, hot rolling, rubber sheets, the performance of lubricants and paints, the aerodynamic extrusion of plastic sheets, metal spinning, and polymer extrusion.Hence, Yashkun et al. 1 evaluated the characteristics of thermal radiation and suction influences on linear shrinking and stretching a sheet under magnetohydrodynamic hybrid nano-fluid for the existence of dual solutions.Waini et al. 2 also discussed the thermal radiation effect on a vertically shrinking sheet with mixed convection magnetohydrodynamic flow in the presence of dust particles and CuAl 2 O 3 nanoparticles.The stretching sheet in the presence of a heat sink/ source, thermal radiation, and chemical reactions of first-order non-Newtonian nanofluid mixed convection flow was addressed by Hayat et al. 3 .Furthermore, Jamaludin et al. 4 reported a numerical investigation of the vertical shrinking/stretching of a sheet under the influence of a heat source/sink, suction, and thermal radiation at the stagnation point in the presence of mixed convection flow and water as the base fluid, and with nanoparticles as titania TiO 2 , copper Cu, and alumina Al 2 O 3 .
Searching for a new nanoparticle that enhances the heat transfer rate with minimum cost is the duty of noway researchers, but it is a challenging task for scientists in the field of fluid dynamics.The suspension of one type of nanoparticles into a nanofluid containing nanoparticles of other metals or non-metals can enhance the thermal properties to obtain a new kind of nanofluid, which is called hybrid nanofluid.In different industrial areas like nuclear system cooling, microelectronics, naval structures, manufacturing, biomedicine, and drug reduction, hybrid nanofluids have many applications.Yashkun et al. 5 used water as the base fluid and copper and aluminum oxide as the nanoparticles, with Joule heating and mixed convection under an exponentially shrinking/stretching surface in the presence of heat transfer.A numerical investigation of unsteady mixed convection, magnetohydrodynamic flow of heat transfer, and hybrid nanofluid pasts of permeable vertical plates is reported by Wahid et al. 6 .Muhammad et al. 7 computed the combined effects of mixed convection, slip, and dissipation due to the heating of hybrid nanofulids formed from Cu and MWCNTs past a surface with a curved-shaped body.Further, many researchers also investigated the effects of hybrid nanofluid flow [8][9][10] .
In varying practical applications such as transportation, electronics, etc., that is, a modern cutting-edge adhesive technologies, such as a flow past stretching a curved sheet have a wide range of applications.Some of the applications are growing crystals structures, plastic sheets preparation, manufacturing of electronic chips and materials, paper industry, cooling process.Due to this case, the curved sheet connected to the flow has come to the attention of more researchers in the present day.Researchers such as Rosca and Pop 11 have reported on an electrical conduction fluid of incompressible, viscous, laminar, two-dimensional, unsteady flow with past shrinking sheets in the presence of a transversely identical magnetic field under a curved surface.Also, Revathi et al. 12 investigated the influence of CH3OH + SiO 2 + Al 2 O 3 , a methanol based hybrid nano-fluid, on a curved stretching surface in the presence of thermal radiation and cross diffusion due to increased activation energy and Brinkman number.The exception to the great output devices of entropy generation, activation energy due to forced natural convection flow over a curved surface of fully developed Darcy-Forchheimer is analyzed by Muhammad et al. 13 .Ahmed et al. 14 developed the phenomena of heat and mass transport under a curving, expanding sheet of Cu-CuO dependent on sodium alginate, NaAlg hybrid nanofluid with boundary layer flow in the presence of nonlinear thermal radiation, chemical reactions, and a magnetic field.Ijaz et al. 15 discussed the behavior of the entropy optimization rate of magnetohydrodynamic flow due to a curved stretching surface in the presence of Joule heating, heat generation, and viscous dissipation.Ibrahim and Gizewu 16 reported the characteristic magnetohydrodynamic flow of a third-order slip boundary condition for a curved stretching sheet under gyrotactic microorganisms and bio-convective entropy generation in the presence of Dufour and Soret effects.
The phenomenon in which energy spreads from a heated surface to its absorption point in all directions in the form of electromagnetic waves is called thermal radiation.The thermal agitation of composite molecules in the body generates thermal radiation.Heating of the room due to an open fireplace is a common example of thermal radiation.The sun, the light bulb, and microwave radiation are typical examples.Furthermore, technological applications of it can be seen in solar power, nuclear power plants, combustion chambers, and chemical processes.The amount of heat energy emitted by the heat surface per unit area is directly proportional to the fourth power of the surface's absolute temperature, as stated by Stefan Boltzmann's law of radiative heat transfer.Wahid et al. 17 evaluated the impacts of radiative mixed convection flow with a hybrid nano-fluid, copper-alumina/water, due to the existence of a magnetic field.Bejawada and Nandeppanavar 18 analyzed the influence of thermal radiation in the presence of micro-polar fluid with magnetohydrodynamic heat transfer flow past a moving porous vertical plate.Lv et al. 19 explored the behavior of activation energy and chemical reactions with two horizontal infinite plates without nano-fluid flow in the presence of variable thermal conductivity, the Hall effect, and thermal radiation past a permeable and stretchable sheet.Hussain et al. 20 analyzed the impact of boundary layer flow on velocity slip, thermal radiation at heat absorption, and magnetohydrodynamics past a permeable exponentially stretching surface.The impact of magnetohydrodynamics, chemical reactions, and thermal radiation, in the presence of melting heat transfer, Brownian motion, and thermophoresis effects is due to a non-linear stretching surface, as discussed by Krishnamurthy et al. 21.For more information, see 22 .
Mixed convection is a blend of free and forced convection in the fluid flow.When the effects of forced flow and free convection/buoyant forces in forced convection are present, a considerable mixed convection flow is produced.Jamaludin et al. 23 presented a theoretical model for hybrid nano-fluids with mixed convection boundary layer flow due to an exponentially shrinking/stretching surface in the presence of injection/suction and viscous dissipation under the porous medium.Gohar et al. 24 explored the influence of mixed convection flow in Casson hybrid nanofluid under a curved past stretching sheet on the contribution, according to Darcy-Forchheimer, of permeable media in an incompressible viscous fluid flow.Ibrahim and Gizewu [25][26][27] examine mixed convection flow with modified Fick's and Fourier's diffusion theories of non-Newtonian tangent hyperbolic fluid in the presence of convective and slip boundary conditions past the non-uniform thickness, bi-directional stretching sheet and entropy generation under a thin film flow boundary condition.
Injection or Suction of fluid through the bounding surfaces, as a consequence, affect the rate of heat transfer from the bounding surfaces, as, for example, in mass transfer cooling, it can significantly change the flow field.In general, injection acts in the opposite manners whereas suction tends to increase the skin-friction and heat transfer coefficients .Injection/withdrawal of fluid through porous cooled or heated surface is of general interest in practical problems involving control of boundary layers, film cooling,etc.This can lead to enhanced cooling /heating of the system and can help from laminar flow to delay the transition [28][29][30] .
Existence of dual solutions in numerical computing have becoming an important topic in fluid dynamics problems for researchers to study its physical significant.Some of the researchers' previous works on dual solutions are:-Shi et al. 31 reported the impacts of the dual solution on Maxwell-fluid flow with mixed convection and thermal radiation in the presence of thermophoresis and Brownian motion past an exponentially shrinking sheet.Yahaya et al. 32 evaluated the characteristics of a dual solution of two-dimensional unsteady mixed convection and thermal radiation flow for a Riga-plate with a stagnation point over convective boundary conditions.Nadeem et al. 33 investigated the behavior of micro-polar hybrid nano-fluids with boundary layer flow past a shrinking/stretching surface by using the Runge-Kutta Butcher method along with the Nachtsheim-Swigert iteration with the existence of dual solutions.Naramgari et al. 34 examined the characteristics of dual solutions in magnetohydrodynamics and mixed convection buoyancy of a non-isothermal stretching sheet in the presence of radiation, injection/suction, chemical reactions, a heat sink/source, and a magnetic field.Mousavi et al. 35 and De et al. 36 discuss the influence of a dual solution on the magnetic hydrodynamics of MgO-Ag/water hybrid nano-fluids of Casson flow past a shrinking/stretching surface in the presence of radiation, suction, and convective boundary conditions.
The aim of this article is to investigate the stability analysis and dual solutions of mixed convection with thermal radiation flow of Cu-Al 2 O 3 /H 2 O hybrid nano-fluid over a curved stretching/shrinking surface due to convective and slip conditions in the presence of suction/injection and magnetic field effects.The findings of the present article will provide fruitful results in the future research in the field of fluid dynamics.To the best of

Problem formulation
In this article, 2D, incompressible, viscous boundary layer flow of hybrid nanofluid with the stability analysis for the existence of dual solutions of mixed convection with thermal radiation flow of Cu-Al 2 O 3 /H 2 O hybrid nano-fluid over a curved stretching/shrinking surface due to convection and slip conditions in the presence of suction/injection are considered.r and s are the coordinates of a radius of curvature R. (u, v) are the velocity components along the s and r directions, and a is the positive constant, where g is the acceleration due to gravity.The magnetic field Bo is taken into account in determining the direction of the flow.A hybrid nano-fluid, in which H 2 O is used as the base fluid and Cu and Al 2 O 3 are used as solid nanoparticles, is taken into considera- tion.Figure 1 presents the physical flow model of the present paper.A system of differential equations under the above assumption for the boundary layers of the flow problem in unsteady form are govern as (see 5,7,11,12,16 ): Equation of continuity: Equations of momentum: Equation of energy: ∂u ∂t www.nature.com/scientificreports/ The boundary conditions defined by 11,16 are as the follows: V w = − √ aν f S , U w = as and > 0, < 0, S > 0, S < 0, 2 and 3 are representing stretching sheet, shrinking sheet, suction, injection, velocity slip and thermal jump parameters respectively.

Steady-state case ( ∂u ∂t = 0)
The subsequent similarity variables employed in the assumed steady-state case of the governing problem in Eqs.
When we compare Eqs. ( 9) and ( 8), we obtain: The dimensionless forms of the associated boundary conditions are: where From Eq. ( 9) P would be written as the following: (5) The skin friction coefficient Cf r and local Nusselt number Nu l (see 37 ) are the physical interests of engineering quantities expressed as follows: where after implementing Eqs. ( 19), ( 15) and ( 16), where are represented by the reduced local Nusselt's number, reduced skin friction coefficient, and local Reynolds number, respectively.
The thermo-physical characteristics of nano-particles and base fluid as given by 5 in the following Table 1: The hybrid nano-fluid under consideration are given by 5 :

Stability analysis
Researchers 38,39 were reported the existence of dual solutions.To determine these solutions are realizable physically, stability analysis of equations ( 11)-( 13) is required.In order to observe which solution is realizable the parameter τ is utilized with an initial value problem where new variables given as the following: ( 14) .
( www.nature.com/scientificreports/One obtains after substitution of ( 19) into the unsteady form of (1) and ( 6) are: subject to: Then, consider the following perturbation functions which is used to check the stability of the steady flow solutions g(ξ ) = g 0 (ξ ) and θ(ξ ) = θ 0 (ξ ) satisfying the boundary value problem ( 11)-( 13) (see 38,39 ) where unknown eigenvalue parameter is γ and the functions G(ξ , τ ) and θ 1 (ξ , τ ) are relatively small compared to g 0 (ξ ) and θ 0 (ξ ) respectively.The stability of the solutions are determined by the sign (positive or negative) of the eigenvalues.Putting of ( 23) into ( 20)- (22) gives the following linearized eigenvalue problem: With their corresponding boundary conditions are: by considering τ = 0 (24)-( 26) are reduced as the following: (20) It should be reported that for various values of Pr the stability of the steady flow solution g 0 (ξ ) is determined by the smallest eigenvalue γ .According to 38,39 , the range of possible eigenvalues can be determined by relaxing a boundary condition on G 0 (ξ ) .Thus, we relax the condition that G ′′ 0 (ξ ) → 0 as ξ → ∞ and for a fixed value of γ ( 27)-( 28) along with the new boundary condition G ′′′ 0 (0) = 1 has to be solved.

Numerical solution
To find the numerical solution of the non-linear ordinary differential equations subject to the boundary conditions of the flow problem, bvp4c solver using Matlab software was employed.To solve the non-linear ordinary differential equations ( 11) and ( 12), with their respective boundary conditions ( 13), we form a system of first order ordinary differential equations by letting in the form of g = f 1 ; Then Eqs.11-13 are reduced to the following form:

Outcomes and discussion
The following values were used to analyze the effects of current governing parameters: M = 0.9, Ec = 1, S = 5, ξ = 0.4, β = −1, K = 50, Rd = 0.2, 1 = 0.1, 2 = 0.2, 3 = 0.2, Pr = 6.2, = −10, φi = 0.05 and φj = 0.02 .The characteristics of stability analysis and dual solutions for the dimensionless parameters of the present study are given as follows: The effect of M; Ec; β ; K; Rd; S; 1 ; 2 ; 3 ; Pr; ; φi and φj on the velocity profile, shear stress profile, temperature profile, skin friction coefficient profile, and local Nusselt number profile are shown in Fig. 2 6 depict the velocity profile g ′ (ξ ) with variations of , S, 1 , and 2 .Figures 2 and 3 illustrate the velocity profile g ′ (ξ ) with variations of and S, respectively.In these figures, boundary layer thick- ness of the first solution is greater than the second solution.As the values of both and S increase, the velocity profile for the first solution decreases and the opposite trend is observed for the second solution.Figure 4 portrays (28)   1 Pr www.nature.com/scientificreports/ the velocity profile g ′ (ξ ) with a variation of 1 .In this figure, the boundary layer thickness of the first solution is smaller than that of second solution.Physically, this is due to the fact that an enhancement in 1 increases the buoyancy force, which results in an increase in the velocity profile for both solutions.Figure 5 labels the velocity profile g ′ (ξ ) with a variation of 2 .In this figure, the boundary layer thickness of the first solution is grater than the second solution.It is physically noted that the velocity profile for first solution is increasing when 2 increases.As a result, the opposite trend was observed for both solutions when various values of 2 were presented.Figure 6  www.nature.com/scientificreports/portrays the velocity profile g ′ (ξ ) with variation of φj.In this figure, the boundary layer thickness of the first solu- tion is greater than the second solution.As the values of φj increase, the velocity profile for the second solution decreases, and the opposite trend is observed for the first solution.Figures 7, 8, 9, 10, 11, 12, 13 and 14 depict the shear stress profile g ′′ (ξ ) with variations of β , K, 1 , 2 , , and S. Figure 7 illustrates the shear stress profile g ′′ (ξ ) with a variation of β .In this figure, the boundary layer thickness of first solution is greater than the second solution.It can be observed that the thickness of the shear stress increases or decreases for the second or first solution.As a result, the shear stress profile increased for the second solution, and the opposite trend was observed for the first solution.Figures 8 and 9 delineate the shear stress profile g ′′ (ξ ) with variations of K and 1 , respectively.It can be observed from both figures that the thick- ness of the shear stress profile for the first solution is larger than that of the second solution, even though after some interval point of ξ the opposite trend was observed.As a result, the shear stress profile increased for the first solution, and the opposite trend was observed for the second solution, though after some interval, the inverse of this statement was observed.Figures 10, 11    stress profile g ′′ (ξ ) with variation of the volume fraction of Al 2 O 3 φj and mass suction parameter S. From Fig. 13 and 14, we came to understand that in both figures, the thickness of the shear stress profile for first solution is smaller than second solution even though after some interval point of ξ the opposite trend was observed.As a result, the shear stress profile decrease for the first solution , and the opposite trend was observed for the second solution, even though some interval points of ξ inverse of this statement were observed.Figures 15, 16, 17, 18, 19, 20, 21, 22 and 23 show the effects of the temperature profile θ(ξ ) with various values of β , Ec, K, , M, φj , Pr, Rd, and S. Figure 15 sketches the temperature profile θ(ξ ) with a variation of β .It is observed from this figure that the thickness of the thermal boundary layer for the second solution is smaller than the first solution.As the values of β increase, the temperature profile for the first solution increases and the opposite trend is observed for the second solution.Figure 16 sketches the temperature profile θ(ξ ) with variation of Ec.It is observed from this figure that the thickness of the thermal boundary layer for the second solution is smaller than the first solution.Physically,the ratio of the dynamic temperature to the temperature or the kinetic energy to the enthalpy driving force for heat transfer is Ec.As a result, as the values of Ec increase, the temperature profiles for both solutions increase.Figure 17 sketches the temperature profile θ(ξ ) with a variation of K.The thickness of the thermal boundary layer becomes thinner as the K impact is increased, even though for the second   solution is larger than the first solution.As the values of K increase the temperature profile of both the first and second solutions decreases.Figure 18 depicts the temperature profile θ(ξ ) with a variation of .It is seen from this figure that the thickness of the thermal boundary layer for the second solution is smaller than the first solution.As a result, The second solution is a decrease; the opposite trend was observed.Figure 19 sketches the temperature profile θ(ξ ) with variations of M. It is seen from this figure that the thermal boundary layer thickness for the first solution is greater than that of the second solution.Physically, the Lorentz force due to the transverse magnetic field has the property of relaxing the fluid velocity and temperature distributions.Accordingly, the temperature distribution of both solutions enlarges as the values of the magnetic parameter M increase.Figure 20 sketches the temperature profile θ(ξ ) with variation of φj .It can be observed that the thermal boundary layer thickness increases or decreases for the second or first solution.As a result, the temperature profile increased for the first solution, and the opposite trend was observed for the second solution via φj .Figure 21 sketches the temperature profile θ(ξ ) with variation of Pr.It is seen from this figure that the thermal boundary layer thickness for the first solution is larger than the second solution.physically, the ratio of momentum diffusivity to thermal diffusivity or kinematic viscosity to thermal diffusivity.From this fact, for various values of the Pr, temperature profile have opposite characteristics for both solutions.Figure 22 delineates the temperature profile θ(ξ ) with variation of Rd.The fact that the relative contribution of conduction heat transfer to thermal radiation transfer is a radiation  parameter.As it can be seen from the figure, the thermal boundary layer thickness for the first solution is larger than that of the second solution.This shows that as the rate of radiation increases, the temperature profile in both solutions also increases.Figure 23 sketches the temperature profile θ(ξ ) with variation of S. It is seen from this figure that the thermal boundary layer thickness for first solution is greater than that of the second solution.As a result, for both solutions, the opposite trend in the suction parameter is increased.
The impact of the dual solution and the critical value for various values of M on the skin friction coefficient and the local Nusselt number are portrayed in Figs.24 and 25 respectively.Figure 24 portrays the behavior of the dual solution and the critical values for M on the skin friction coefficient.Critical value for M = 12.0, 12.1, and 12.2 are c1 = −1.141414,c2 = −1.292929, and c3 = respectively, where < 0 up to which the solution exists.Dual solution exist between ci < where i = 1,2 and 3.For the value of M uphill, of the absolute value of ci increase for both the skin friction coefficient and the local Nusselt number.The same property, which can be easily understood is registered in Fig. 25, the characteristic of dual solution and critical value for various value of M on the local Nusselt number.3 and 4 show the effect of dual solutions with varying M, S, Br, and Pr values on the g ′′ (0) , Cfs, and Nus of Cu − Al 2 O 3 /H 2 O , a hybrid nano-fluid.The results present the properties of the non-dimensional param- eters M, S, Br, and Pr on g ′′ (0) , Cfs, and Nus in an easy-to-understand manner, as shown in Tables 3 and 4. For different estimations of Pr the smallest eigenvalue are shown in Table 5.From this Table 5 we generalize that the first solution is stable (physically realizable) whereas the second solution is unstable (not physically realizable).

Conclusions
In this study, two-dimensional, incompressible, viscous boundary layer flow of hybrid nanofluid with stability analysis and a dual solution of mixed convection and thermal radiation with hybrid nano-fluid over curved surfaces past a stretching/shrinking surface are considered under slip conditions.A hybrid nanofluid,in which water is used as the base fluid and copper and alumina are used as nano-particles, and magnetic field are taken into account.The bvp4c method is used to solve the numerical solution of the governing equations of the present study.The main outputs of this study are: www.nature.com/scientificreports/ • For the values of both and S, the velocity profile for the first solution decreased, and the opposite trend was observed for the second solution.
• As the values of φj increase, the velocity profile for the second solution decreases, and the opposite trend is observed for the first solution.
• As the values of β increase, the shear stress profile increases for the second solution, and the opposite trend is observed for the first solution.
• For the values of K and 1 , the shear stress profile increased for the first solution, and the opposite trend was observed for the second solution, though after some interval points, the of this statement was observed.www.nature.com/scientificreports/ • For the values of K, the temperature profiles increase, and for both the first and second solutions, they decrease.
• As the values of Ec increase, the temperature profiles for both solutions increase.
• For the values of S, Pr, Rd, M, and , the upwind thermal boundary layer of the first solution is larger than the second solution.• For the value of M uphill, the estimation of the absolute value of ci increases for both the skin friction coef- ficient and the local Nusselt number.• In the second solution, increasing the values of β , Pr, S, Ec, and M has a similar effect on g ′′ (0) , Cfs, -θ ′ (0) ,and Nus.• In the first solution, increasing the values of S, and Pr on g ′′ (0) and Cfs result in a decrease.
• The first solution has a positive eigenvalue, whereas the second solution has a negative eigenvalue.

Figure 1 .
Figure 1.Diagram of the present work.
and 12 label the shear stress profile g ′′ (ξ ) with variations of 2, and φi .It can be observed that from all Figs. 10, 11 and 12 the shear stress profile increased for the first solution at some interval point of ξ , and after this interval, inverse of this statement was observed.Even though the shear stress profile decreases for the second solution in all Figs. 10, 11 and 12. Figures 13 and 14 illustrate the shear

Figure 25 .Table 5 .
Figure 25.Portrays the existence of dual solutions on Re −1/2 s Nu s with variation of M.

Table 1 .
Thermophysical characteristics of nanoparticles and base fluid.

Table 3 .
portrays the effect of the dual solutions, which are the first and second branch solutions of M, Ec, S, Pr, and β respectively, on g ′′ (0) and Cfs.

Table 4 .
portrays the effect of the dual solutions, which are the first and second branch solutions of M, Ec, S, Pr, and β respectively, on -θ ′ (0) and Nus.

Table 2
summarizes the current result and compares it with previous work in the open literature for various values of K on -Cfs.The result confirms good agreement, this shows that the present method gives an acceptable result.Tables